SMOOTH APPROXIMATIONS FOR FRACTIONAL AND 
MULTIFRACTIONAL FIELDS 
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Abstract. We construct absolute continuous stochastic processes that con- 
verge to anisotropic fractional and multifractional Brownian sheets in Besov- 
type spaces. 



1. Introduction 

Fractional Brownian motion (fBm) with Hurst parameter H is a continuous 
centered Gaussian process with the covariance function 

EB,B, = ^{t^" + s'"-\t-s\'"), t,s>0. 

It has stationary increments that exhibit a property of long-range dependence for 
H > 1/2, which makes fBm a popular and efficient model for long-range dependent 
processes in Internet traffic, stock markets, etc. 

For numerous applications one needs a multi-parameter model for a long range 
dependence. For most of those applications, including image processing, geophysics, 
oceanography, etc, two parameters are enough. 

There are different possibilities to define a two-parameter fractional Brownian 
motion, or fractional Brownian sheet (IBs). One is so-called isotropic fractional 
Brownian sheet with the covariance function 

EBtBs = i {\\tf" + pf ^ - \\t - sf") , t,se Rl. 

As the name suggests, its properties are the same in all directions, which is not the 
case in many applications, especially in image processing. For such applications a 
better model is anisotropic fractional Brownian sheet with the covariance function 

EB,Bs = 1 n + - 1^. - s.\'"^) , t,se Rl. 

i=l,2 

It also has stationary increments in the sense that the distribution of (Bt+s — Bt) 
does not depend on t. 

But the stationarity of increments of fractional Brownian process and sheet 
means that the behavior of them is the same at each point, and this substan- 
tially restricts the area of their application. In particular, they does not allow one 
to model situations, where the regularity at a point depends on the point, as well 
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as the range of dependence. In view of this, recently several multifractional gen- 
eralizations of fBm and fBs were proposed in order to overcome these limitations 
such as moving average multifractional Brownian motion (mBm) [B^ , harmonisable 
mBm j2j , Volterra-type mBm [8j , Different multiparameter extensions of mBm were 
studied in [Tl[3l|4]. 

In this paper, we work in a general setting by considering a Gaussian random 
field B{t) in the plane which is continuous almost surely and satisfies the following 
condition on its increments: for all s,t G [OiTi] x [0,T2] 



where C > and A > 1 are some constants. This, in particular, includes anisotropic 
fractional and multifractional Brownian sheets. 

Our main goal is to construct approximations for B(t) in a certain Besov, or 
fractional Sobolev, spaces, by absolutely continuous fields. This allows one to ap- 
proximate stochastic integrals with respect to fractional Brownian sheet by usual 
integrals and consequently, to approximate solutions of stochastic partial differen- 
tial equations involving fractional noise by solutions of partial differential equations 
with a random source, which in many aspects are similar to non-random partial 
differential equations. 

The paper is organized as follows. In Section [2] we give necessary definitions 
concerning Besov spaces. In Section [3] the main result. Theorem 13. li about the 
convergence of absolutely continuous approximations of the field B(t) in Besov space 
in probability, is proved. In Section |4] we study fractional and multifractional 
Brownian sheets. We show in particular that these fields satisfy the conditions of 
Theorem 13. II and can be approximated by absolutely continuous random fields. 



In this section, we define functional spaces which are similar to spaces of Holder 
continuous functions. They play an important role in definition and analysis of 
stochastic integrals with respect to fractional random fields (see e. g. [5J Lemma 



Let s,t e M^, s = (si,S2), t = {ti,t2)- Write s < t ii si < ti and S2 < h- 
For s < t denote [s,t] = [si,ti] x [52,^2] C M^. Let T ^ {Ti,T2) e (0,00)^, 
[0,r] = [0,Ti] X [0,r2]. For a function /: R we consider two-parameter 

increments 



Asf{t) := f{t) - /(si, t2) - /(ii, S2) + /(s), s, t e R\. 
Consider a function /: [0,r] M. For /3 e (0, 1)^ and t e [0,T] we denote 



E(B(S1, S2) - B{si,t2) - B{ti,S2) + BihM)? < Ci\h - Sil \t2 - S2|)\ 



2. Definitions and notation 



2.2.16]). 





\f{t)\ + vi^im + (^^^(/)(t) + ^i^"^{f){t). 
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Let W^^'^^^ = Wo^'''^' ([0, T]) be a spaee of measurable functions / : [0, T] R, with 

ll/lloA,/3. = sup ^f^^W<oo, 
*e[o,T] 

Define also wf''^' = Wf'''''([0,T]) as a space of measurable functions /: [0,T] 
R, with 



1 /•*M/t-(«,S2)-/t_(s)| 



|AJ(t)| 



1 



dv 



+ , , JV'-'I < CO. 



/ 



M (^i-«i)'+^^(^2-S2)i+'5^ y 
where /t_(s) := f{s) - /(si,t2-) - /(ii-,S2) + /(*-)• 

3. Main result 

Let {Bt,t e [0,T]} be a random field which satisfies the following conditions 

(1) Bt is a Gaussian field; 

(2) there exists constants C > and A > 1 such that for all s,t G [0, T] 
(1) £{^sBtf <C{\h- s^\\t2- S2\f; 

(3) the trajectories of Bt are continuous with probability one. 

One example of such field is anisotropic fractional Brownian sheet (see Introduc- 
tion). 

We consider the approximation 

^t = - / B,ds = - / B,+tds. 



Theorem 3.1. For all Pi, ^2 G (0, A/2) 
Proof. Denote 

AB| := B| - = 1 / {B^+t - Bt) du. 

Then 

K{ABl) = AB| - ABf^,,^ - AB|^,,^ + AB| = 1 / (B„+t - 

— Bui+Sl,U2+t2 + Bsi,t2 ~ Bui+ti,U2 + S2 + Bti, S2 + — iJg) rfU. 

According to the Cauchy-Schwarz inequality, 

(2) E(As(Ai?j ))^ < ^ / E (-B„+t — Bt — Bu.^+si,u2+t2 + -^si.ts 
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Considering dU), we obtain 

E(A,(ABf))2 < 4 / {E{Au+sBu+tf + E(A,BO') du 
(3) ' 



< 4 / 2C{\h-Si\\t2-S2\)Uu^ACi\h~.Si\\t2-S2\ 
[O.e 



On the other hand, ^ imphes 

'[0,£ 



E(A,(Ai?n)' < ^ / (E(A 

£ J[0,£]2 



(4) /■ / 



8C 



|ti-si|V|i2-s2ne^ 



A + 1 

because, in view of ([1]), 

E(A„,+,,,t,B„+()2 <C\h- sil^ u^, 

E(A„^+Sj_S2i3„j+t^^„2+S2)^ < C \tl — Sl\ U2, 

E{A,B^,,+,,^t2f <Cu^\t2-S2\\ 

Let (5 G (0, A — 2max{/3i,/32}). We study three cases. 
Case 1: |ti — si| < e. Based on the estimate ([3]), we get 

E(A,(AS|))2 < 4C(|ti - sil |i2 - 52|)^ < 4CT|(|ti - si| |i2 - S2|)^"'e'. 
Case i^; |ti — si| > e, |t2 — S2I < e. Similarly to the Case 1, we obtain 

E(A,(AB|))2 < 4C(|ti - sil |i2 - S2\)^ < 4CT/(|ii - si| 1^2 - S2|)^~'e'. 
Case 3: \ti — si| > e, 1^2 — S2I > £■ Based on the estimate (|4]), we get 

E(A,(ASnf < 3^ (|ti - + 1*2 - .2I') 

< 4C (r/ + r|) - sii |i2 - s2|)^-'^£^ 

Thus, in all 3 cases we have 

E(A,(AS[))2 < 4C (r/ + r|) - sil \t2 - s2\)^-'e'. 
Since As{ABf) has a normal distribution, then for p > we have 

(5) E |A,(ABnr < Ciilh - 5i| 1*2 - S2|)(^~^)^/"£'f/^ 

where Ci = Ci{p,Ti,T2) = 

According to the two- parameter Garsia-Rodemich-Rumsey inequality ([71 The- 
orem 2.1]), for all p > 0, Qfi > p~^, Oi2 > there exists a constant C2 = 
C2(ai,a2,p) > such that 

(6) \As{ABt)f < C2 iti - sir^p-i \t2 - S2r^p-^ c, 



SMOOTH APPROXIMATIONS FOR FRACTIONAL AND MULTIFRACTIONAL FIELDS 5 



where 

|A.(ABn|P 



[o,T]2 \xi - yiT'^^^ - y2\ 



— dx dy. 



We choose < 6l < (A - 2max{/3i, (32} - S)/2, p = |, ai = as = Then, 
taking into account (O, we obtain 



E|A.(Ai?g)|^ 
[o,T]2 \xi - yiT'^'*'^ \x2 - ?/2r 



J[0,T]2 



= Cie^P/^ [ dxdy^CiT^T^e 

J[0,T]2 

Therefore, from ([6]) we get 



where C3 = CiC2T^Tl 

The last estimate imphes that for any k G (0,1) there exists such that prob- 
abihty of the event 



A, := {for all s,t e [0,T] : |A.(Ai3f)| < c„(|ti - 5i| \t2 - S2|)(^-29-5)/2^V2| 



is not less than 1 — k. 

As we have chosen < (5 < A — 2 max {/3i, /32} and < < — max {/3i, /32}, 
then - /?! - 6* > 0, /i2 := ^ - /32 - 6* > 0. 

By the definition of the norm ||-||-^ we have 



sup 



|A,(Ai3f; 



0<s<r<T \{tl~S^Y'it2~S2f^ 

, 1 r*^ |A.(Ai?g^,j| 1 Z-*-^ |A.(Ai3r„J| 

(<2 - S2Y- A, - + (ti - 4 (t; - 52)1+'^^ 

, / |A.(Ai?.^)| 
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Therefore at the set 



lAB" 



< sup c^e 

0<s<t<T \ 



<5/2 



\'l2-l 



dv 



5/2 



sup 

0<s<t<T 



h2-l 



dr 



1 1 

hi /i2 



hih-2 



Then for any a > 



hm P ( HAS 



> a I < K, 



because for sufficiently smaU s one has Cne^/'^Ti^T2^{l 
Hence, when k 0+ we have 



e-i-0+ 



(« 



Km P HAS'' 



> a U= 0. 



□ 



Remark 3.2. The convergence in probabiUty in the last theorem may be not suf- 
ficient for some applications. For example, one may want to use this theorem to 
get an approximate solution to stochastic PDE with a fractional noise by solving 
a usual PDE with random force. She is not able of course, to solve for all w's and 
takes some fixed uj. She knows, of course, that there is a subsequence of solutions 
converging almost surely, but apriori it is not known which subsequence is it. So 
another subsequence (depending, say, on lo) may converge to something different or 
there might be several such subsequences. To overcome this problem, we give below 
a proof that for e„ = 2~" one has an almost sure convergence, and give moreover 
an estimate for the rate of convergence. 



Let 



a{e) 
bis) 



£5p/2 1^1+7 £^ ^ > 0. 



|A.(ABf)| 



([7]) implies that 
Then 



a[en) 



V a(2-") 
^ 6(2-") 

n>l ^ ' 



< 



— ^ 0, n ^ oo. 



n>l 
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Therefore, 

a(e„) 



— i> 0, n — i> oo, a. s. 



Hence, 

a(£„) < C(w)6(e„) a.s. 
So for all s, t G [0, T] and n > 1 

|A,(AB-)| < Ci/^'(c.)4/2ln(i+^)/^'e„(|ti - .i| \t, - a.s. 

Using the definition of the norm ||-||-^ we obtain 

l|A5fJli,....<^^/W^^T^ (l + ^+zJ^ + ^i 

X e„ ^ 0, n ^ C50, a. s. 

4. Examples 

4.1. Fractional Brownian sheet. Fractional Brownian fields in the plane can 
be defined in various ways. We consider the so-called anisotropic random fields 
that possess the fractional Brownian property coordinate-wise (see e.g. [5j Sec- 
tion 1.20]). 

Definition 4.1. A random field \Bf G [0,r]} is called a fractional Brownian 
sheet with Hurst index H = (ifi, H2) e (0, 1)^ if 

(1) Bp isa. Gaussian field such that Bp ^Q,te dM.\] 



(2) EBf = 0, EBf Bf = 1 n (tr- + - 1^. 

1=1,2 ^ 



\2H, 



This field has a continuous modification. Its increments satisfy the equality 

E(A,i3f)' = |ti-si|'^^ \t2-s2r'. 

Hence, for B^ the inequality ([T]) holds with A = 2min{i/i, _ff2}. Therefore, ac- 
cording to Theorem 13. 1[ for Hi g (^,1) for aU /3i,/32 G (0,-ffi A H2) one has a 
convergence of approximations 

\\b"^' ~B"\l Ao, s-^0+, 

II IIJ-,Pl,P2 ' 



BP'^XI I B"ds. 



where 



4.2. Mult ifr actional Brownian sheet. We consider a function 

H{t) = (i?i(i),i?2(<)): [0,T] ^ (1/2,1)2. 
Let fj,, V be constants such that 

- < u < min i/i(t) < max Hi(t) <v< 1. 

2 te[o,T] te[o,T] 

Assume that there exist positive constants ci, C2 such that for all i, s e [0,r] 
(HI) |i7,(t) - ff,(s)| < ci (|ii - siT + 1*2 - S2r), 

(H2) \^sH,{t)\<C2{\ti~Si\\t2~S2\Y. 
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Definition 4.2. Multifractional Brownian sheet ^B^'^*\t e [Oi^jj with Hurst 
function H{t) is defined as 



dWu, 



^H,(t)-l/2" 

' i=l,2 

t e [0, T], where s+ — max {s, 0}, W — {Ws,s G M^} is a standard Wiener field. 



Denote Yt = Bf^*\ 

Theorem 4.3. The trajectories of are continuous with probability one. 

H (t^ 

Proof. We prove that trajectories ofYt — B^ are continuous with probabihty one 
on any rectangle [a,b] C (0,r] such that ||5 — a\\ < 6. Moreover, according to 
Lemma 2.2], the constant 5 > can be chosen such small that for all s,t G [a,b] 
the following inequality holds 



E (Yt - Y,f <CiJ2 \ti - Si I'" < 2Ci \\t - 



As Yt is a Gaussian field, then for any a > there exists C2 > such that 

Taking a > 2/ we obtain that the field Yt is continuous on [a, 6] with probability 
one according to Kolmogorov theorem. □ 

Theorem 4.4. There exists a constant C > such that for all s, t G [0, T] 

(8) £{AsYtf <C{\t^-s^\\t2-S2\f'' 
Proof. Denote s' ~ {si,t2), t' = (ii,S2). Then 

AsYt ^Yt- Yt, + Y, - Y,, =Ai+A2 + A3 + A4, 

where 

A2 = B^^'^ - Bf (*) + Sf (^') - B^^''\ 
A3 = i?f (*) - + - (^'^ 

A4 = B^^'^ - - Sf (*) + Bf (*'). 

Hence, 

(9) E(A,rO' < 4(EA? + Ev42 + EAj + EAj). 

We estimate each of 4 terms. 

As the random field B^ is a fractional Brownian motion when H{t) = const, 
then 

(10) EAl < C3 1^1 - \t, - s,f"'('^ < Ci{\ti - si\ \t, - s.lf^. 
Consider A2. 

A2 = A21 +A22, 
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where 



i-ui) 



{t2 - U2) 



H'At)- 



H2{t)- 



(si - 

[Si - Ui)^ 



(si - Ml) 



Hi(s')- 



Hi{s')-i 



du 



H2(t)-^ 



( [(t2 - «2)f - (-^^2)f ^*^"^] - [(S2 - U2) 

By Lemma [5TH 

Ev42^ < X2(i?l(t) - i/i(.s'))' • Ki \t2 ~ S2|'^^(*) 
Taking into account Condition (E(T]), wc get 

(11) ^Al^<C^{\h-Si\\t2-S2\f^ 

Now we estimate ^22 • 



i-U2) 



H2{t)-h 



(si - Ui) 
([(^2 -W2) 



+ 



(S2-U2) + 

"/, -.li2(s')-l 
[t2 - U2) + 



' - (-U2 
{-U2f/'^-'' 

( ^H2(s')- 



dui 



{S2 - M2) + 



dU2 



( ^H2(s' 
{-U2) + 



By Lemma [ 

^Al^ < Kisl"'^''^K^{t2 - S2f^{H2{t) - H2{s')f 

Considering the Condition (Etl]), we obtain 

(12) EAl^<C^{\tl-Si\\t2-S2\f^' 

(HH and (HI]) imply that 

(13) EAl < C^{\tl ~ si\\t2 ~ S2\f'' 
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Further, by Lemma 15.41 

EAl < L [{Hi{t) - Hi{t') + Hi{s) - Hi{s')f + (i/aW - i?2(<') 
+ H^i-s) ~ H2{s')f + {{H,{t) - H,{t')f + {H2{t) - H2{t')f 
+ {H,{s) H,{s')f + {H^is) HM)f) {{Hiii) H,{s')f 
+ [H-2{t) - H2{s')f + {H,{s) - H,{t')f + {H2{s) - H2{t'))')] , 
and taking into account the conditions ([T]) and ([2|), we get 

(14) EAl<Csi\ti-si\\t2-S2\f^. 
The estimation of A4 is analogous to that of A2. We have 

(15) EAl<Cc>{\h~si\\t2~S2\f^. 

Combining ^, we get (HI). □ 

Theorems 14.31 and 14.41 imply that multifractional Brownian sheet i?^^*-* satisfies 
the conditions of Theorem 13. II Therefore, 



Corollary 4.5. For any /3i,/32 G (0,/^) 



H{t),e 



H{t) 



where 



H{t),e 



l,/3i,/3: 
ti+e pt2+e 



ds. 



5. Appendix 

In this section we prove some technical lemmas that have been used in the proof 
of Theorem 14.41 



5.1. Bounds for fractional Brownian motion. Let (i7,J^, P) be a complete 
probability space. 

Definition 5.1. Fractional Brownian motion with Hurst index H £ (0,1) is a 
centered Gaussian process — {Z^,t > 0} with stationary increments and the 
covariance function 

1 



\t- 



i2H\ 



It is well-known that fractional Brownian motion has a continuous modification 
and can be represented in the following form. 



Z" - Ch 



(t-<-^-(-<- 



where W is a Wiener process, Ch — '^^rui+i%f'^'^ (^^^ Chapter 1.3]) 



Let i < /!< i/min < ffn 



:'{H+l/2) 

< v < I. Consider a family of random variabiles 



^H- 



dWu 



ie[o,r],He[Ai,H. 

Lemma 5.2. There exist positive constants Ki, K2, K3 such that 
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(1) for all tiM e [0,T], H G [i/,„i„, i?,„ax] 



(16) 



,2H 



(2) for all t e [0,T], i/i,ff2 e [i^min, i?max] 

(17) E (zf^ - Zf^)' < i^2(i/i - H2f. 

(3) /or an h,t2 e [0,T], i/i,i/2 £ [ffmin,ffmax] 

(18) E (z/f^ - z^^^ - z^- + ^) ' < i^3(ti - ^2)'^(i^l - i^2)' 

Proof, (i) By the definition, 



which entails ([TB]) because Cj^^ is bounded when iJ G [/i, v\. 

(ii) The inequality (fTT]) is a corollary of (fTSl) . 

(iii) We prove p^ . Let for definiteness, t2 < ti, H2 < Hi. We can write 



"(tl-^^)f ^-(t2-^i)f ^-(tl-u) + '"^+(t2-w) + ' 
= /l+/2, 



where 
h = 

h = 



ih - u)"^-i - {t2 - U)"'-^ - {h - uf'-^ + {t2 - U) 



du. 



't2 



du. 



By the theorem on finite increments, there exists h G [iJ2,^^i] such that 

Ji = {Hi - H2Y f ' \{ti - uf-^ \n{ti -u)- {t2 - u)''--^ \n{t2 - u) 

(19) ^"^7 

^{Hi-H2f j\v-uf-"^ (^+[h-^^\n{v-u)^dv 



du 



First we prove that for ahy h E[pL,v\ 
rt2 r rti 



(20) 



{v ~ uf-'^/^dv 



du<C^{ti-t2f>', 
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where C3 > is a constant. Actually 



t2 



-I 2 



du 



< 



r2t2—ti 


f\v- 


uf-^'^dv 


' — 00 


VJt2 




r2t2-ti 




- uf~^'^dv 


J —00 


Jt2 





du 



du 



2t2-tl 
t-2 



{v - uf-'^/'^dv 



du 



H 2 



du 



^{h~t2)"'{i2~2h)-' + {h-l/2)-^) 

^ (Itt) + " ^'^^"^ + " '/^^"'^ 

< (ii - i2)'^(r + 1) ((2 - 2^.)-! + (/i - l/2)-2) . 



It is obvious that for all e, (5 > there exist positive constants and such 
that 



(21) 
(22) 



|lna;| < C^x^' 
\liix\ < Csx^ 



if < x < 1, 
if X > 1. 



First we show that 



(23) 



t2-l 



{v - m)''-3/2 {l + {h- 1/2) ln(i; - u)) dv 



-I 2 



du < C^i-t^ff. 



Applying the inequality ([22|) with 5 := {v — iJ,„ax)/2, we obtain 



{v ~ uf-'i {l + {h- 1/2) ln(u - u)) dv 



du 



< 2 



t2-l 



*2 



r.t2-l 



du 



2{v - l/2YCi 



ti 



{v - uf-i+^dv 



t2 



du. 



Using we get 

Secondly we prove that 



(24) 



t2 



t2-l 



{v - uf-^/^ {l + ih~ 1/2) ln(w - u)) dv 



du<Cz{ti-t2f''. 



Choose £ := (/i — i/mi„)/2. Then (PT|) implies that if u G \t2 — 1, ^2], v G [<2, ^i] then 



|ln(w — u)\ 



In(r + 1) + In 



T+ 1 



< ln(r + 1) + Ce 



T+ 1 
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Therefore 

t2 
t2-l 



tl 



{v ~ m)''-3/2 (1 + (/i - 1/2) ln{v - u)) dv 

rt2 



du 



< 2(l + (!.-l/2)ln(T+l))-^ 



t2-l 



du 



2{v-l/2YCt{T + 1) 



\2e 



t2-l 



[v - uf-^/'^-'dv 



t2 



du. 



Considering (|20l), we obtain (p4)) . 
Combining ([T9l). (pS]). pi]), we get 



It remains to estimate l2- By the theorem on finite increments, there exists 
h e [i?2, Hi] such that 



/2 = {Hi - H^f / \{h - uf-^'^ \niti - u) 



du. 



Choose £ :— — iJi„in)/2. Then (pi]) implies that if u e [i2,ii] then 



which entails 

h < 2(-ffi - H^ 



ln(r + I) + In 



r + 1 



< in(r + l) + Ce 



tl — u 

T+1 



<C7(il-t2)'n^l-^2)'. 

Now the proof is complete. 

5.2. Bounds for integrals. Let i < /i < i?niin < ^^max < < 1- 
Lemma 5.3. Let 

f{t,U, h) = {t- uf-^'^ - {-u)X-^^\ te[Q,T],ue^,he [i?,„i„,-ffmax]. 

Then for all t <E [0,T], he [i?min, -ffmax] 



□ 



(25) 
(26) 
(27) 



[fit,u, h)fdu < +00, 
[fhh{tiU,h)fdu < +0O. 
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Proof. We prove (|27|) . Inequalities psp and ([26l) are proved in a similar way. 







{v — u) 







- i ) ln^(w - u) + 2 ln(w - u) ) 



du 



+ / [t - uf'-Hn^t - u)du =: h + h + h, 



where 



{v - uf^"^ ( [h- - ] \n{v - u) + 21n(v - u)] dv 



h = 



[v - uf~^ [ [h--] \n{v - u) + 21n(w -u)]dv 



du. 



du, 



h = [ {t - u)^''-Hn\t - u)du. 
Jo 

We study each of three terms. 

1. Applying the inequality (P^. we obtain 



h < 



f {v-uf -^/^ (^(^h-^^C]iv~u)^^ + 2Cs{v-u)^^dv du 



<2 h 



CI 



du 



8Cl 



du < Ct^^" < oo 



(the last estimate follows from the inequality ([20]) ). 
2. Consider /2. 



h = 



T+ 1 



r + 1 



h--] ln^(r+ 1) + 21n(T+ 1) ) dv 



du. 



Using inequality (PT|) . we get 



h < 



2e 



+ 2Ce{T + lf{v - u)-' + 1) 



+ 21n(r+ 1) ) dv 
where the last estimate follows from the inequality (|20|) . 



n 2 



du < Ct^^' < oo, 
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3. Consider I3. 



h= / v^'''Hn'' vdv < / v^''-Hn^vdv+ / v^''~Hn\dv 
Jo Jo Jl 

< I v^i"-^ hi' vdv + / v^'^-^ln'vdv < 00. 







Thus, inequality ([27]) holds. □ 

5.3. Bounds for fractional Brownian sheet. Let ^ < fi < i/min < ^max < v < 
1. We consider a family of random variables 

Bt / ( (<i -ui)+ - (-wi)+ 



X ((t2 - ^i2)f - {-U2)+-^'"') dWu, 

t e [0, T], H, H' e[fi,v],i^ 1, 2, where W = {W,, s e R^} is a Wiener field. 

Lemma 5.4. There exists a constant L > such that for all t G [0,r], Hi,H[ G 
[^^min, ^^max]? * = 1, 2, 3, 4, the following inequality holds 

(28) E [b"' "'' - Bf + b"' "'^ - Bf'^^'^y 

< L {{Hi -H2 + H3- Hif + {H[ -H'^ + H'^ - H'^f 
+ {{H, ~ H^f + {H[ - H'^f + (i/3 - H^f + m - H',f) 

X ((i/i - + {H[ - H'^f + {H2 - + {H'2 - H'^f)) . 

Proof. Denote 

f{t,u,h) = {t-u)X-"^ -{-ut-^'\ 

.fi(h) = f{ti,ui,h), f2{h) = f{t2,U2,h). 

Then 

E (Bf'^"'' - b"^'"'' + b"''"'' - Bf''"''Y = [ F^{t, u) du, 

F{t,u) = f,{H,)f2{H[) - fi{H2)f2{H'^) + fi{H3)f2{H',) - f,{H,)f2{H',). 

Consider two cases. 

Case 1: {Hi ~ i?2)(i?3 - H^) > 0. 
In this case 

\Hi-H2\ <\Hi-H2 + H3~Hi\, 
\H3-Hi\ < \Hi-H2 + H3~Hi\. 

We have 

F{t, u) ^Fi{t,u)+ F2{t, u) + F3{t., u), 

where 

Fi{t,u) = {fi{Hi) - fi{H2))f2{H^), 

F2{t,u) = (/i(i/3) - fl{Hi))f2{H'^), 

F3{t,u) = fi{Hi){f2{H[) - f2{H'^) + fi{H,){f2{H'^) - f2{Hi)). 
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By the mean value theorem, there exist hi G [Hi A H2,Hi V H2] and /12 G [-^3 A 
Hi, Hz y Hi] such that 

h{Hi)- h{H2)= f[{hi){Hi~H2), 
fi{Hs)- h{Hi)^ f[ih2){Hs~Hi). 

Applying Lemma 15.31 we get 

/ Fi{t,u)du<{Hi^H2)' f {f[{hi)fdu[ {h{H'^)fdu, 
< Ci{Hi - H2 + H:i - Hif . 
In much the same way, we have 

/ Fi {t, u)du < C2iHi - H2 + Hi - Hif . 

JR2 

Hence, 

/ F'^{t,u)du<3( Ff{t,u)du+ Fi{t,u)du+ F^{t,u)du] 

7r2 \Jr2 Jr2 Jr2 J 

< CaiHi -H2+H3- Hif + 3 / u) du. 

JR2 

Consider the latter term. There are two possible cases. 
Case la: {H[ - H^)iH^ - H'^) > 0. 
In this case 

\H'i - H2\ < \H[ - H2 + H!^ - H'^ \ , 
\H^ - H'^l < \H[ - H2 + H'^ - H'^\ . 

By the mean value theorem, there exist /13 £ [H[ A H2,H[ V H2] and /14 e [H^ A 
H'i.H'^y H'^] such that 

f2{H[)-f2{H'2)^mM-H'2), 
h{H'3)-f2{H'i)=r2{hi){H'3-H'i). 

Using Lemma 15.31 we estimate 

/ Fi{t,u)du<{H[-H'2f [ {h{Hi)fdu[ {mz)fdu 

JR2 JR2 Jr2 

+ {H'3 - H'^f [ {h{Hi)f du [ (mi))' du 

JR2 JR2 

<Ci{H[-H!2+H!,-H',f. 
Case lb: {H[ - - H'^) < 0. 

Without loss of generality, we assume that — iJjl > \H^ — H'^\ and H[ < H2 
(hence, H'i< H^). 

Put Hi ^ H'2 + H'i - Hl^. It is not hard to see that Hi e [H[,H'2]. 



Fs{t,u) ^ F3i{t,u) + F32{t,u), 
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such that 



where 

F3lit,u) = (/2(i/0 - /2 (i/l)) , 

F32{t,u) = (/2 (i/i) - /2(i?^)) + h{Hi){h{H'^) - fM))- 

By the mean value theorem, there exists /15 G , Hi 

F3i{t,u) = /i(i/i)/^(/i5) (i/i - ffi) = h{Hi)f^{h^){H[ -H^+H!,- H',). 
By Lemma [Ol 

/ < (F; -H^+H',- Hif [ {h{Hi)f du [ {r^{h,)f du 

We estimate ^32(^,11). 

F32(i, u) = - ^ ' ' /i(-ffi)/^ (i/i + dx 

+ / ' ' h{Hi)f^{H'^ + x)dx 

pH'^-H'^ /■Hi 

< / / fi{y)f2[Hi+x] dydx 



fi{H4)f2{x + y) dydx. 

Hi 

By the mean value theorem, there exist xi, 2:2 £ [0, -^3— -^4], J/i £ [H1AH4, Hi\/H4], 
2/2 e Hi A H!^,HiW H'^ such that 

fHs-H'i /■Hi 

= (ifi - il4)(i/^ - H'4)f'i{yi)f2 {Hi + , 
/ f'i{Hi)f^{x + y)dydx 

'Hi 

= (ffl - Hi) {H'^ - Hi)f[iH4)f^{x2 + 2/2) 
= {H!, - H'2){H'^ - H'4)f[{H,)f^'{x2 + y2). 
Using Lemma 15.31 we get 

/ i^322(^, n)du < C^{H'^ - H'^fiHi - H^f + Cj{H'^ - H'^fiH'^ - H'^f . 



Thus, in the case lb (|28l) holds. 
Case 2: {Hi - iJ2)(i?3 - H^) < 0. 

Without loss of generality, we assume that \Hi — H2\ > 1^3 — i?4| and Hi < H2 
(hence, H4 < H3). 

Put Hi ^ H2 + H4 - H3. It is not hard to see that Hi e [Hi,H2]. 
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We have 

F{t,u) = Gi{t,u) + G2{t,u) + G3{t,u), 

where 

G2{t,u) = - h{H2))h{H[) + iMHs) - h{Hi))f2iHi), 

G3{t,u) = MH2){f2{H[) - f2{H^)) + /l(F3)(/2(if^) - f2{Hi)). 

Terms Gi{t,u) and G2{t,u) are estimated similarly to Fsi{t,u) and Fs2{t,u) in 
Case lb. We have 

[ Glit, u)du <CsiHi-H2+H3-Hif, 

/ Glit, u)du < Ca{H3 - Hif{H[ - H'^f + C^iH^ - H^fiH^ - H2f. 

It remains to estimate G3{t,u). We consider two cases. 

Case 2a: {H[ - H'^){H'^ ~ H^) > 0. 

G3{t,u) can be estimated similarly to F3{t,u) in Case la. We get 
/ Glit, u)du < CniH[ -H^+H^- H'^f. 

Case 2b: - iJ2)(^^3 - H^) < 0. 

This case can bo considered in a similar way to Case lb. Without loss of gener- 
ality, we assume that — iljl > 1-^3 ~ ^4 1 ^^'^ H[ < Then we obtain 

/ Glit, u)du < Ci2{H[ -H'^ + H'^- H'^f 
Jr2 

+ C,3iH'3 - H'^fiH3 - H2f + Ci4(i?^ - H'^fiH'3 - H!,f. 
Thus, now the proof is complete. □ 
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